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Abstract 

To describe the scattering of superconducting quasiparticles from non- 
magnetic (Zn) or magnetic (Ni) impurities in optimally doped high T c 
cuprates, we propose an effective Anderson model Hamiltonian of a localized 
electron hybridizing with d x 2_ y 2-wme BCS type superconducting quasiparti- 
cles with an attractive scalar potential at the impurity site. Due to the strong 
local antiferromagnetic couplings between the original Cu ions and their near- 
est neighbors, the localized electron in the Ni-doped materials is assumed to 
be on the impurity sites, while in the Zn-doped materials the localized elec- 
tron is distributed over the four nearest neighbor sites of the impurities with 
a dominant d x 2_ y 2 symmetric form of the wave function. Since both scatter- 
ings from the localized electron and the scalar potential are relevant, localized 
resonant states due to their interplay are formed below the maximal super- 
conducting gap. With Ni impurities, two resonant states are formed above 
the Fermi level in the local density of states at the impurity site, while for 
Zn impurities a sharp resonant peak below the Fermi level dominates in the 
local density of states at the Zn site, accompanied by a small and broad reso- 
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nant state above the Fermi level mainly induced by the potential scattering. 
This is exactly what has been observed in the scanning tunneling microscopy 
experiments. In both cases, there are no Kondo screening effects. From the 
calculated spin relaxation functions, we find that the 3d localized electron in 
both Ni and Zn doped materials displays a weak magnetic oscillation. This 
result is consistent with the signal of a spin- 1/2 magnetic moment exhibited 
by nuclear magnetic resonance measurements in YBa2Cu306+<5 doped with 
Zn or Ni impurities. The local density of states and their spatial distribution 
at the dominant resonant energy around the substituted impurities are calcu- 
lated for both cases, and they are in good agreement with the experimental 
results of scanning tunneling microscopy in Bi2Sr2CaCu208+«5 with Zn or Ni 
impurities, respectively. Thus the scanning tunneling and nuclear magnetic 
resonance experiments on Zn and Ni substituted cuprates are interpreted self- 
consistently in a unified fashion. 

PACS numbers: 71.10.Hf, 71.27.+a, 71.55.-i, 75.20.Hr 
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I. INTRODUCTION 



Recent experiments on <i-wave superconducting cuprates have shown many interesting 
phenomena in the presence of magnetic (Ni) and non-magnetic (Zn) impurities substituting 
the Cu ions in the copper-oxygen plane. A series of nuclear magnetic resonance (NMR) 
experiments on optimally doped YBa 2 Cu30 6+( 5 with impurities have shown that each 
non-magnetic impurity (Zn) induces a local spin- 1/2 magnetic moment on the neighboring 
Cu sites, while a magnetic (Ni) impurity gives rise to a localized spin- 1/2 magnetic moment 
on the impurity site. On the other hand, scanning tunneling microscopy (STM) has been 
used to probe the quasiparticle scattering around single Zn or Ni impurity ions in optimally 
doped Bi 2 Sr 2 CaCu208+<5 with high spatial and energy resolution |H|. In the obtained STM 
spectra, a very sharp localized resonance peak is found just below the zero-bias voltage 
on the Zn impurity sites |4j], while two localized resonant states are formed well above the 
zero-bias voltage on the Ni impurity sites ||. The spatial dependence of the local density 
of states (LDOS) in the vicinity of the individual impurity ions reveals the characteristic 
features of d x 2_ y 2 wave superconducting (dSC) states, but the pattern is different in these 
two cases. 

A large number of theoretical papers have been devoted to the interpretation of these 
beautiful experimental data. According to one school of thoughts, this is due to a potential 
scattering, either spin-independent, or magnetic (classical magnetic moment) jBHlOf. In fact, 
the quasi-bound state around a potential scatter in c?-wave superconductors was discussed 
before the experiments were done ||||]. On the other hand, there is an alternative explana- 
tion ascribing the resonant scattering as due to the Kondo effect in ci-wave superconductors 



1113, extending earlier studies of the Kondo effect in systems with reduced density of 



states at the Fermi level fll4yi5| . It is fair, however, to say that several important issues re- 
main unexplained, for example, the strong particle-hole asymmetry, the difference between 
the Zn- and Ni-doped cases, concerning both spectral distribution and spatial pattern. Thus 
a full understanding of these experiments is still lacking. 
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In an earlier communication [[T(J, we proposed an effective Anderson Hamiltonian with 
hybridization of a localized electron with quasiparticles in c?-wave superconducting states 
to describe the scattering on quantum magnetic impurities. In the limit of the Hubbard 
repulsion U —>■ oo, we have shown the existence of a sharp resonance above the Fermi 
level and absence of the Kondo screening. Moreover, a "marginal Fermi liquid" behavior 
was found for the impurity electron self-energy in the strong coupling limit. In this paper 
we extend our study to include the potential scattering effect on the impurity site. The 
interplay of the quantum magnetic scattering and potential scattering turns out to play a 
very important role in determining the spectral distribution and spatial pattern. 



The NMR and neutron scattering |T7j experiments show that Zn doping has a 



stronger disturbance to the antiferromagnetic background in cuprates than Ni doping. There 
is a kind of common understanding that Zn doping induces a local magnetic moment dis- 
tributed mostly over the nearest neighbor Cu ions, whereas Ni impurity has an underscreened 
S — 1/2 moment, sitting on the Ni ion itself. This picture is consistent with the electron con- 
figurations. Zn 2+ ions have a closed shell (3c? 10 , S = 0), with an extra c?-electron compared 
with Cu 2+ which is not entirely sitting on the Zn site. The electronic structure studies |Tj| 
do show that the electron charge density is higher on neighboring oxygen sites and nearest 
neighbor Cu 2+ ions. On the other hand, Ni 2+ ions (3c? 8 , S = 1) have an extra hole of spin 1/2 
on site, compared with Cu 2+ ions, and the deficit of charge is concentrated on the Ni site. [18 



In this paper, we propose a unified model to describe these two cases. We consider a local 
charge carrier (3c? electron in the Zn case and 3c? hole in the Ni case, to be called Anderson 
electrons for short, and such a treatment is allowed due to the particle-hole symmetry of the 
BCS state) with a strong on-site Coulomb repulsion U, is hybridizing with superconducting 
quasiparticles of d x 2_ y 2 symmetry. In the Ni-doped case, the Anderson electron is sitting 
on the impurity site itself, whereas in the Zn-doped case, it is distributed over the nearest 
neighbor Cu 2+ ions. The dominant contribution comes from the linear combination with 
d x 2_ y 2 symmetry. Moreover, we assume an attractive 5-type potential scattering in both 



cases, which can be justified to some extent within the three-band Hubbard model . JLS 



In the above effective model, due to a strong interplay between the quantum impurity 
scattering (described by the Anderson model) and the potential scattering of the quasipar- 
ticles in dSC state, we find a sharp resonance state below the Fermi energy mainly due 
to the Anderson electron, and a broad potential scattering peak in the positive energy in 
the LDOS for the Zn-doped case. In contrast, two resonance peaks due to Anderson elec- 
tron and potential scattering are formed well above the Fermi level for the Ni-doped case. 
These distinct features fully agree with STM experiments @,|]. Moreover, the spatial pat- 
tern of the resonance states differs in these two cases, also in agreement with experiments. 
The calculated dynamical spin response functions show typical structures corresponding to 
the existence of local magnetic moments. These results demonstrate convincingly that the 
proposed unified model describes the essential physics of these experiments self-consistently. 

The rest of the paper is organized as follows. In section II, the nickel impurity substitution 
is considered, with both Anderson electron and potential scattering center located at the 
impurity site itself. In section III, we discuss the interplay of quantum magnetic scattering 
on the Anderson electron distributed over nearest neighbors and on-site potential scattering, 
corresponding to the Zn-doped case. Finally, concluding remarks are made in section IV. 

II. MODEL CALCULATION FOR NICKEL IMPURITY SUBSTITUTION 

To construct a simple model, we assume that a BCS-type weak coupling theory is ap- 
plicable as a phenomenological framework for high-T c optimally doped superconductors, 
though the underlying mechanisms are very different. When the Cu ions on the CuC>2 plane 
are replaced by magnetic Ni impurities, there will be one extra 3d hole on each impurity ion. 
Due to the particle-hole symmetry of the dSC quasiparticles, we can work in the electron 
representation, and assume that the magnetic Ni impurity is effectively described by the 
Anderson localized electron with a strong Hubbard repulsion. Moreover, a scalar attractive 
potential at the impurity site has been added to the model Hamiltonian as caused by in- 
teractions with the nearest neighbor Cu ions. When the correlations between the magnetic 
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impurities on different sites are ignored, the model Hamiltonian is defined as 



u = e k4 CT c kff + J2 A k (ckt c -ki + h - c - 
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where e k = h 2 k 2 /(2m) — tp is the dispersion of the normal state electrons, A k = A cos2# 
is the <iSC order parameter with the gap amplitude Ao, and W < is the strength of the 
attractive potential. When Nambu spinors are introduced 
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we rewrite the model Hamiltonian in a matrix from 
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where a z and a x are Pauli matrices. The equations of motion for the dSC quasiparticles 
yield the generalized T-matrix exactly 



f(iu n ) = V 2 [a z - WG^iu^Ga (iu n ) [a z - WGo^)]- 1 
+W[a z - WG {iu n )\-\ 



(3) 



where G (iu n ) = jr J2k — e kO"z — A k cr z ) 1 is the GF of the dSC quasiparticles at the 
impurity site. When the scalar potential is absent, the T-matrix is reduced to 



T (iu n ) = V a z G d (iu n ) a z , 



(4) 



which was obtained previously in the single impurity Anderson model in dSC electron state 
|T6|. On the other hand, if only the scalar potential is present, it becomes 
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f (iu n ) = [(Waz)- 1 - d Q (iu n )]-\ (5) 

which is the same result as derived by Balatsky et al. . From the generalized T-matrix, one 
can clearly see that both hybridization with the localized electron and potential scattering are 
equally important in determining the low energy excitations. At zero temperature, analytical 
continuation is used to calculate the perturbed GF through the following relation: 

G (r, r'; u) = 6° (r - r', u) + G° (r, u) f (u) G° (-r', u) . 

The LDOS of the dSC quasiparticles around the impurity is given by 

N (r, u) = Im Gn (r, r; uo) . 

We would emphasize that the above relation between the GF of the dSC quasiparticles and 
the Anderson electron is exact, and the only quantity to determine is Gd(u), which will be 
calculated approximately. 

When we take the infinite U limit, the localized electron operator is expressed as (p + = 
f+~b, fib + J i n the slave-boson representation [2D|,ZT|], where the fermion f a and the boson 



b describe the singly occupied and hole states, respectively. The constraint b + b+Yl fa fa — 1 
has to be imposed. When a mean field approximation is applied, the boson operators b and 
b + are replaced by a c- number bo, and the constraint is treated by introducing a chemical 
potential Ao- Therefore, the mean field Hamiltonian is written as 

k ^ V k,k' 

+ ~7Tf (^k °z4> + h.c.^j + e d <j) + a z (j) 



N k 

+e d + X (b 2 - 1). (6) 

where <p + = ( J+ ; ^denotes the Nambu spinors of the localized electron and the renor- 
malized parameters q = e^ + Xo an d V = b$V . It has been known that the slave-boson mean 
field approximation produces the correct low energy physics of the conventional Anderson 
impurity model in the entire parameter range at zero temperature, though strong quantum 



fluctuations are present near the Fermi surface from the gapless excitations of the conduc- 
tion electrons ||22j| . In the present model, we are dealing with dSC quasiparticles, where 
a superconducting gap exists on most part of the Fermi surface except four nodal points, 
we thus expect that the slave-boson mean field approximation should work better than the 
conventional case, giving rise to more reliable results at zero temperature. 

The standard techniques lead to Gf(iuj n 
energy of the localized electron becomes, 



iu n - l d o z - £/ (iu. 



where the self- 



E f (tuj n ) = V 2 [G^(iu n ) - Wa z }-\ 
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At T = 0, the ground-state energy changes due to the presence of localized electron and the 
scalar potential are: 
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where D is the band width of the normal conduction electrons, a(ui) = 
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. The corresponding saddle-point equations are derived as 

2u 2 f3(l + b 2 (3) - 2{W/V 2 )u 2 af3 [l d - b 2 (W /V 2 )uj 2 a(3] 
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n cu 2 (1 + %PY + [e d - b 2 (W/V 2 )u 2 af3] 
For given parameters D, A , T = ttNfV 2 , e d , and c = 7rATp|W^|, one can obtain saddle 
point values of b and Aq. In the following calculations, we choose A as the energy unit, 
D/A = 20, and T/A = 0.2. 

The phase diagram can be obtained for a specific value of the potential scattering strength 
c, similar to the phase diagram given previously flUfl . When e d is less than a threshold value, 
6o is zero, leading to a free local magnetic moment decoupled from the dSC quasiparti- 
cles. Above the threshold value of e d , b^ rises sharply and then saturates. The low energy 
behavior of the model is controlled by the nearly empty orbital state, characterized by a 
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strong hybridization between the Anderson electron and the dSC quasiparticles. In the hole 
representation, the empty orbital phase corresponds to a single hole occupation at the Ni 
impurity site. Such a situation can not be realized in the conventional Anderson impurity 
model, where the low energy behavior is always controlled by the Kondo screening fixed point 
in the single electron occupied phase. The dSC quasiparticles actually provide us with such 
a possibility to explore new low-energy physics in the empty orbital regime — the strong 
coupling behavior in the single hole occupied regime. In a sense, such a strong coupling 
behavior is also universal, corresponding to a new strong coupling fixed point without any 
characteristic energy scale. However, the scalar potential will introduce some extra features 
to the LDOS around the impurities. 

We have seen that both scalar potential and localized electron scatterings are relevant in 
the dSC state, and two resonant states are formed with a dominant sharp localized resonance 
above the Fermi level. In Fig. 1, we plot the DOS for the localized electron (u) with 
different values of for two potential scattering strengths c = 1 and 2. For a small e^, the 
resonant peak positions can be roughly determined as 



Clearly, the former resonance peak is mainly induced by the localized Anderson electron, 
while the latter weaker resonance is caused by the attractive potential. As moves up, the 
magnitude of the first localized resonance f^ es decreases and becomes broader with shifting 
upward of the peak position, while the magnitude of the second resonance Q^es emerges 
as a pronounced peak at a relative high energy. When ~ 0.30, the profile of the two 
consecutive resonances agrees well with the measured DOS at individual Ni site ||. This 
does also confirm our assumption of the attractive scalar potential at the Ni impurity site, 
because the corresponding localized resonant state would appear below the Fermi energy if 
the scalar potential scattering is repulsive. 

The nature of the double resonant peaks in the DOS of localized electrons can also be 
understood from the imaginary part of its retarded self-energy: 




(10) 
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ImZ)' 1 L + iQ + ) =Im- J£ , (11) 

which has been plotted in Fig. 2a. In contrast to the case in the absence of attractive 
scalar potential (solid line), a single peak appears above the Fermi energy (dashed line), 
corresponding to the pole of Eq. ([[!]). This peak becomes narrower with increasing 
the potential strength c (dash-dotted line). Moreover, the existence of the scalar poten- 
tial greatly suppresses the spectrum weight of the self-energy far away from the Fermi 



level. However, the marginal Fermi liquid behavior |23| obtained previously is preserved: 
ImEj; 1 H ~ -b 2 T Q£)as u -> 0. 

Furthermore, in order to reveal the characteristic behavior in the regime of the empty 
orbital (single hole) localized state, we calculate the dynamic spin response function of the 
localized electron 

XdM = ((Si® : S" d {t))) (12) 

corresponding to particle-hole collective excitations, from which the spin relaxation function 
Sd(u) = ~~ Xd(a ^'° ' can be evaluated. Sd(u>) is related to the neutron scattering cross- 
section and the calculated value is displayed in Fig. 2b. There is a pair of symmetric peaks 
demonstrating a weak magnetic oscillation for the 3d localized hole, in contrast to the uni- 
versal result of one central resonant peak obtained in the Kondo regime of the conventional 



Anderson impurity model [p4}| . Such a behavior is consistent with a local spin 1/2 magnetic 
moment on the Ni impurity site observed in the NMR measurements 0. We also notice 
that closer the localized resonant state is to the Fermi level, stronger magnetic oscillation the 
localized hole exhibits. This can be understood by considering the linear energy dependence 
of the LDOS for the dSC quasiparticles: the localized resonant state becomes more and more 
delocalized as one moves away from the Fermi level. As for the self-energy, the presence 
of potential scattering (c ^ 0) reduces substantially the magnitude of the spin relaxation 
function. 

Now consider the LDOS of the dSC quasiparticles N (r, u) in the spatial range < r < £. 
In high T c dSC state, £ is about 10A, or roughly 3 lattice spacings. In Fig. 3, the LDOS vs 
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frequency is shown for r=0.07£ from the impurity site along the vertical (or horizontal) and 
diagonal directions in real space. Similar to the LDOS for the localized electron, there are 
two resonances above the Fermi energy along the diagonal (the gap nodal) directions and 
the LDOS is entirely hole-like. In contrast, for vertical and horizontal (the gap maximum) 
directions, two more resonant peaks appear below the Fermi level as well, which are slightly 
asymmetric in the line shape and strongly reminiscent of the localized bound states induced 



by paramagnetic impurities in the conventional s-wave superconducting state ||25|| . In both 
of these cases, increasing the potential strength will sharpen and enhance the localized reso- 
nances significantly. As the impurity energy level e<j moves up, the quasiparticle resonances 
become broader, exhibiting a similar dependence as the resonance of the Anderson electron 
itself. 

It is very important to calculate the spatial distribution of the LDOS of dSC quasipar- 
ticles to compare with corresponding STM results. The LDOS around the impurity at the 
first resonance energy is displayed in Fig. 4a for the hole-like excitations to = — Ql es and in 
Fig. 4b for the electron-like excitations u = Q} Tes as a function of spatial variables for c = 1,2, 
and 3 with e^/Ao = —0.05 in a logarithmic intensity scale. The quasiparticle resonances 
induced by substituted impurities are highly localized around the impurity, and the spatial 
oscillation of these resonant states is visible. The largest amplitude of the quasiparticle 
resonance occurs in the neighborhood of the impurity, and the local electronic structures 
distinctly differ in Fig. 4a from that in Fig. 4b. For to = Ql es , the LDOS exhibits a four- 
fold symmetry along the directions of the gap nodes for all distances, consistent with the 
dSC quasiparticles. For u = — f2* es , the LDOS is strongly enhanced in the gap maxima 
directions at distances r <C £. Further away from the impurity (r~£), it is confined to the 
neighborhood of the diagonal directions, preserving the four-fold symmetry. At the same 
time, the LDOS of the dSC quasiparticles at the impurity site (the center of the diagram) 
and around the impurity site get enhanced by increasing the strength of the scalar potential. 
In fact, these features have been seen in the STM experiment 0, namely, the resonance at 
positive energy is enhanced in the nodal directions, while the resonance at negative energy 
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is more visible in the anti-nodal directions. To the best of our knowledge, these distinctive 
features have not been explained by other theoretical models considered so far. The absence 
of particle-hole symmetry is natural in our consideration, due to the quantum origin of the 
resonant states, and it is in full agreement with STM experiments. |J 

III. MODEL CALCULATION FOR ZINC IMPURITY SUBSTITUTION 

In the high T c cuprates, the substitution of Cu ions by Zn impurities adds an additional 
3d electron to each Cu ion. By carefully analyzing the experimental results [l],@,|J, one 
can notice that local antiferromagnetic correlations around the non-magnetic impurities get 
enhanced. As a result, there should be a localized magnetic moment distributed mainly 
over the four nearest neighbors of Zn impurities with an attractive scalar potential at the 
impurity. In the following, we will work in the electron representation, and assume that 
the localized state is still described by the Anderson localized electron with a strong on-site 
Hubbard repulsion. The effective model Hamiltonian is thus given by: 

V.a k 

W + V- ,+ i 

iV k,k',cr cr 

j v * (ctJa + h.c.) + Ud^d^d^ (13) 



N k(T 

Compared with the previous model Hamiltonian for the nickel substitution, the momentum 
dependence of the hybridization strength takes the form of Vk = l / cos26 l , because the 
Anderson electron sitting on the four nearest neighbor sites of the Zn impurity hybridizes 
with the dSC quasiparticles with a dominant d x 2_ y 2 symmetric form of the wave function. 



This is the most effective coupling channel in the dSC [|TT],|T3[. The model Hamiltonian now 
describes an extended Anderson quantum "impurity" coupled with the dSC quasiparticles. 
In terms of the Nambu spinors, the GF of the conduction electrons is deduced from the 
equations of motion 
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G(k,k';iu n ) = 5 K ^G (k,iuj n ) + — G (k, iw n )T k y(iw n )G (k', iuj n ), (14) 

where the generalized T-matrix becomes momentum dependent. After Fourier transforma- 
tion, the above equation is expressed explicitly as 

G (r, r'; iu n ) 

= G° (r - r', iu n ) + G° (r, iu n ) % (iu n ) G° (-r', iu n ) 
+V 2 G° t (r, iu n ) a z G d {iu n ) cr z G° t (-r', iu n ) 
+V 2 G° (r, iw n ) a z G d (iu n ) a z G° (r , iu n ) G° p (iu n ) G° (-r', iw n ) 
+\/ 2 G° (r, iw n ) G° p (iu n ) G° (t , iu n ) a z G d (iu n ) (-r', iw n ) , (15) 

with a momentum independent part of the T-matrix 

f (icu n ) = G° p (icu n ) + V 2 G° p (iu n ) G° (r , iu n ) a z G d (iu n ) (r z G° (r , iw n ) G° p (iu n ) , 

where G p (iu n ) = W[a z — WGo(iuu n )Y l is the resulting GF for the o?SC quasiparti- 
cles at the impurity site in the presence of the potential scattering only, G° (r ,iu n ) = 
j^v VkG (k : iu n ) corresponds to the GF at the localized electron site, and G° (r,iu n ) = 
J2k VkGo(k, iu n )e lk ' r is the corresponding Fourier transform of the quasiparticle GF with 
a form factor of the localized electron. In the present continuum model, the positions of the 
nearest neighbor of the impurity r will be set according to our chosen parameters. 

When the slave-boson mean field approximation is applied , the quasiparticle GF for the 
localized electron is derived as Gf(icu n ) = 
ing self-energy is given by 

t f (iw n ) = V 2 a z [G° s (r 0: iu; n ) + G° (r ,iu n ) G° p (icu n ) G° (r , ^ n )]"V,, (16) 

where G° s (r ,iiu n ) = ^Ek^o(k,«w n ). In order to carry out further calculations, the 
various GFs introduced above are simplified as follows 

G (iu n ) = -iu n {-jr-^j a(«^n), 
13 
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. v/l + ^/Ao) 2 ; yJl + (u/A Y V V 1 + (^/ A o 
where and -E'(x) are the complete elliptical integrals of the first and second kind, 

respectively. At zero temperature, the ground-state energy changes due to the presence of 
Zn impurity and scalar potential are: 

D 



5e = J du In ^< 
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l + ^HAH + e d -b 2 A 2 (cu 



where D is the band width of the normal conduction electrons, and 



A 1 (u) 
A 2 {u) 



L + (cAp) 2 a(uj)j{uj) 
T 2 + c 2 u 2 a 2 (uj) ' 
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The saddle-point equations are thus derived as 
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(21) 
(22) 



(J ' uj 2 [1 + bfo^A^u)] 2 + [e d - blA 2 {u)f 

Notations here are the same as used in Eqs.(8) and (9). These equations are solved 
self-consistently to obtain the phase diagram which is essentially the same as for the Ni- 
doped case, i.e. the physics is dominated by the single-hole occupied state. Again, due to 
the presence of both attractive scalar potential and quantum impurity scatterings, localized 
resonant states in the quasiparticle spectrum are induced. However, unlike the Ni-doped 
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case, a dominant sharp resonance associated with the localized Anderson electron is formed 
above the Fermi level, while a small and broad resonance state associated with the scalar 
potential appears below the Fermi level. This difference is mainly due to the different 
location of the Anderson electron: on-site vs nearest neighbors. We plot the DOS for the 
localized electron Nj (u) with different values of for two potential scattering strengths 
c = 2 and 5 in Fig.5. As moves up, the magnitude of the resonance above the Fermi level 
slightly decreases and the peak position is shifted away from the Fermi energy. 

To compare with results for Ni substitution, the imaginary part of the retarded self- 
energy for the localized electron is calculated for the strengths of the potential scattering 
c = 1,2,5 and e^/Ao = —0.2, displayed in Fig. 6a. The main difference from the Ni doping 
case is that a single narrow resonance peak shows up below the Fermi level and becomes 
sharper as the strength of potential scattering is increasing. In Fig. 6b, the spin relaxation 
function for the localized electron Sd(u)) is also shown. A magnetic oscillation is displayed 
and the amplitude is one order of magnitude larger than that in the Ni case (Please note the 
solid line corresponds to c = 1 now) , implying the single 3c? hole on the nearest neighbors of 
the Zn impurity atom does display a substantial local magnetic moment! Certainly, we 
do not claim here a full explanation for the local magnetic moment formation induced by 
Zn impurity doping. However, our calculations do indicate the existence of a local magnetic 
moment on the nearest neighbors of the Zn impurity atom. 

To confront the STM experimental results, we calculate the LDOS of the dSC quasipar- 
ticles around the Zn impurity sites under e^/Ao = —0.3,-0.1 for two strengths of scalar 
potential c = 2,4, displayed in Fig. 7. The LDOS on the Zn impurity site has clearly shown 
that there are a very sharp resonance just below the Fermi level and a small and broad 
resonance above the Fermi level. The former is induced by the hybridization between the 
localized hole with the dSC quasiparticles, while the latter is caused by the potential scat- 
tering. Moreover, in contrast to the Ni substitution result, the DOS near the gap edges are 
completely suppressed to zero, which is consistent with the fact that the superconductivity 
may be suppressed by Zn impurities stronger than that by the Ni impurities. On the nearest 
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neighbor sites r = (±a, 0) and (0, ±a), where a = 0.05£ and £ = Hvf/Aq is the coherence 
length of dSC state, two resonant peaks appear below and above the Fermi level, respec- 
tively. However, on the next nearest neighbor sites ro = ±(a, a), the LDOS obtained is very 
much similar to that on the Zn impurity site. All these features are in good agreement with 
the STM experimental results [Q. 

Finally, the spatial distribution of the LDOS of the dSC quasiparticles is calculated. The 
LDOS around the Zn impurity at the main resonance energy is displayed in Fig. 8a for the 
hole-like excitations u = —Vt res and in Fig.8b for uo = Q res for the electron-like excitations 
as a function of spatial variables for c = 1,2, and 3 with e^/Ao = —0.05 in a logarithmic 
intensity scale. The quasiparticle resonances induced by the doped non-magnetic impurity 
are extended from the impurity to one coherence length of dSC. The largest amplitude of the 
quasiparticle resonance occurs at the neighborhood of the impurity, particularly at the next 
nearest neighbors, and the local electronic structures at distances r <C £ drastically differ in 
Fig. 8a from Fig. 8b. However, further away from the impurity site, both structures of the 
LDOS are almost the same, where the enhanced LDOS is extended along the gap maxima 
directions. Moreover, the spatial LDOS of the dSC quasiparticles is not enhanced as the 
strength of the scalar potential is increasing. The LDOS pattern in the left column for c = 1 
case is closer to the distribution of LDOS obtained from the STM spectra [|J. These results 
seem to show that the strength of potential scattering in the present model Hamiltonian, 
needed to describe the experimental data is not as strong as assumed in previous studies 
with purely potential scattering @,§]. 

IV. CONCLUSIONS 

In this paper, we propose an effective model Hamiltonian of an Anderson localized elec- 
tron (hole) hybridizing with d x 2_ y 2-wave BCS type superconducting quasiparticles supple- 
mented with an attractive scalar potential at the impurity site. Within a unified framework, 
we can describe the non-magnetic (Zn) and magnetic (Ni) impurity scatterings for cuprate 



16 



superconducting quasiparticles of optimally doped high T c materials. 

Due to the nature of (i-wave superconducting quasiparticles, both quantum scattering on 
localized electron and scalar potential scatterings are relevant interactions in the low energy 
limit, showing a strong interplay with each other and leading to localized resonant states 
below the maximal superconducting gap. In the Ni case, two localized resonant states are 
formed above the Fermi level in the LDOS at the doped impurity site, while in the Zn case 
a sharp resonant peak caused by the localized hole is found just below the Fermi level to 
dominate the LDOS at the Zn impurity site, accompanied by a small and broad bound state 
above the Fermi level induced by the potential scattering. These characteristic features are 
in exact agreement with experiments ||||. Moreover, the spatial patterns of the resonant 
states in Zn and Ni doped cases are quite different in STM experiments. It is significant 
that these differences can be reproduced easily in our simple model. The key assumption 
is that the 3d electron is located on the Ni site, while it is distributed over nearest neighbor 
Cu ions in the Zn-doped case due to the local antiferromagnetic background. The excellent 
agreement with experiments convincingly demonstrates the correctness of this assumption 
and the self-consistency of our model. We would emphasize these are generic features of our 
model, not due to adjustable parameters in the theory to fit the experiments. 

We should point out that in both cases, there are no Kondo screening effects on the 
local magnetic moment due to the localized electrons, i.e., the local magnetic moments do 
not couple to the superconducting quasiparticles in the low energy limit, and the dominant 
behavior is the single hole occupied state, characterized by a strong hybridization between 
the localized and delocalized electrons. Moreover, from the calculations of the spin relaxation 
functions for the localized electron in both Ni and Zn doped materials, we find that the 
additional 3d localized hole displays a weak magnetic oscillation. Such a magnetic behavior 
is consistent with the spin-1/2 magnetic moment observed in NMR measurements. 

We would like to mention that our results do not contradict those of the Kondo model 
study for systems with reduced density of states at the Fermi level |]IT|- rr3|jr5|J2B|f27j ] or d- 



wave superconductors. The important factor in those studies is the presence of a strong 
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particle-hole asymmetry, due to, for example, the potential scattering. The behavior of 
systems controlled by the "strong asymmetric coupling" fixed point is very similar to what 
we have obtained here. In particular, there is no single Kondo resonance peak at the Fermi 
level in that case, either. Instead, the resonance is shifted from the Fermi level by the 
energy of the order Tr. In our opinion, identifying this resonance peak to the standard 
Kondo screening effect is somewhat misleading. The NMR experiments in Ni, Zn-doped 
cuprates samples ||l|j2lj28H can be interpreted using the "Kondo temperature" Tk as a char- 
acteristic energy scale. Recent experiment seems to show such a scale exists even in 
the superconducting state. On the contrary, the SQUID measurements do not exhibit any 
signature of such a "Kondo" scale. || We believe the resistivity measurements would be able 
to identify contributions from the Kondo scattering. However, our paper and related stud- 
ies pT|- |13| , p!5[p6| , p7| demonstrate that in low temperature superconducting state one should 



not expect the standard Fermi liquid behavior, corresponding to the Kondo strong-coupling 
limit. To our feeling, the effective Anderson model used in our paper has the advantage of 
treating different physical situations in a unified self-consistent fashion. 

A few related remarks are in order: 
(1) In this paper we did not address the issue of the magnetic moment formation in non- 



magnetic doping case. A lot of efforts have been devoted to this problem ||29|| , but it is still 
not fully understood, in particular, the issue of spinon confinement which is the key argu- 
ment for the "proof, is entirely open. What we have been doing is to assume the existence 
of the local moment, described by the Anderson model and find out that the final results 
are consistent with the starting assumption. 

(2) To reconcile the discrepancy of the potential scattering model || with the STM experi- 
ments on Zn-doped BSCO, namely, the presence of a strong resonance peak at the Zn site 
in experiments, the argument of tunneling matrix element between layers (so-called "filter" 
effect) was suggested. We have not included this effect in our calculations to confront 

the STM experiments. To our knowledge, this issue is still open, and the Kondo model 



calculation 11.13] seems to show the same — irrelevance of this effect. 



(3) We have not considered the effect of disorder on the quasiparticle density of states in 
dSC which may not vanish at the Fermi level. This is still an issue under debate. |5T| Our 
hope is that the main conclusions of our paper will survive in the limit of dilute doping. 
The case of strong disorder is beyond the scope of our present work. 

(4) The interplay of magnetic and non-magnetic scatterings is always important for super- 
conducting states even in the s-wave case p5 |. It is true that the potential scattering alone 
cannot give rise to a bound state inside the gap in that case (Anderson theorem). How- 
ever, the presence of the potential scattering will affect the position of the bound state due 
to the magnetic moment scattering (see first paper of Ref. [p5||). What we have shown is 
that the scattering behavior in the <i-wave nodal direction is quite different from the s-wave 
superconductor, whereas in the gap maximum direction it is rather similar. For the low- 
energy physics, the quantum magnetic scattering is the dominant effect. Finally, we would 
emphasize once again that our result is valid only in the strong coupling, large ?7-limit of 
the Anderson impurity in a ci-wave superconductor, which is very different from that for the 
conventional normal metallic state. 
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Figures Captions 



Fig. 1. The DOS N d (to) of the localized electron, (a) for c = 1 and (b) for c = 2, with 
e d /A = —0.05, 0.05, 0.15, and 0.30, denoted by solid, dashed, dotted, and dash-dotted 
lines, respectively. 

Fig. 2. The imaginary part of the self-energy for the localized electron is delineated in 
(a), and the corresponding local spin relaxation function S d (ui) in (b), for e d /A = —0.05 
and for c = 0, 1 and 2, represented by solid, dashed, and dash-dotted lines, respectively. 

Fig. 3. The LDOS iV c (r=0.07£, u) of the dSC quasiparticles for different potential 
scattering strengths c = 2 and 5, with e d / A = 0.0 (solid line), and 0.3 (dashed lines) in 
units of Np. (a) and (c) are along the directions of the gap maxima, while (b) and (d) along 
the directions of the gap nodes. 

Fig. 4. The spatial distributions of the conduction electron LDOS around the Ni impurity 
at the first resonant energy for c = 1,2 and 3. (a) u = —^l es and (b) u> = +^l es - Here 
e d /A = —0.05, and a logarithmic intensity scale is used. The ranges of x and y in each 
subplot are set from — 0.6£ to +0.6£, and £ is the coherence length. 

Fig. 5. The DOS of the localized electron N d (cu). (a) corresponds to c = 2, and (b) to 
c = 5, with e d /A = —0.4, —0.2, and 0.0, denoted by solid, dashed, and dash-dotted lines, 
respectively. 

Fig. 6. The imaginary part of the self-energy for the localized electron is displayed in 
(a), and the corresponding spin relaxation function S d (u) in (b) for e d /A = —0.2 and for 
c = 1,2 and 5, denoted by solid, dashed, and dash-dotted lines, respectively. 

Fig. 7. The LDOS of the dSC quasiparticles on the Zn impurity site, on the nearest 
neighbor sites, and on the next nearest neighbor sites with e d /A = —0.3 (solid line), and 
—0.1 (dashed lines) in units of N F . (a) corresponds to c = 2, and (b) to c = 4. 

Fig. 8. The spatial distributions of the LDOS of dSC quasiparticles around the Zn 
impurity at the resonant energy for c = 1,2 and 3. (a) uj = —Q res and (b) u = +fi res - 
Here e d /A = —0.2 and a logarithmic intensity scale is used. The ranges of x and y in each 
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subplot are set from —0.3 to +0.3 in units of the coherent length £. 
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